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We study the dynamics of priority-queue networks, generalizations of the binary interacting pri- 
ority queue model introduced by Oliveira and Vazquez [Physica A 388, 187 (2009)]. We found 
that the original AND-type protocol for interacting tasks is not scalable for the queue networks 
with loops because the dynamics becomes frozen due to the priority conflicts. We then consider 
a scalable interaction protocol, an OR-type one, and examine the effects of the network topology 
and the number of queues on the waiting time distributions of the priority-queue networks, finding 
that they exhibit power-law tails in all cases considered, yet with model-dependent power-law expo- 
nents. We also show that the synchronicity in task executions, giving rise to priority conflicts in the 
priority-queue networks, is a relevant factor in the queue dynamics that can change the power-law 
exponent of the waiting time distribution. 
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In the last century, queueing theory has proved useful 
for various problems ranging from operations research 
to telecommunications There is a recent resurgence 
of interest for the queueing theory among the statisti- 
cal physics community with the application to the prob- 
lems in human dynamics. Specifically, various queue- 
ing models based on the prioritization of tasks, or the 
priority queue models to be short, have been intro- 
duced to account for the heavy-tailed distributions ob- 
served in the waiting time and response time distribu- 
tions [1, H, 0, H, [E 0) Hi- The priority queue model is 
grounded on the assumption that the human dynamics 
is the result of an inherent decision-making process of 
the individual, with implicit priorities assigned for every 
tasks in his/her task queue, according to which he/she 
decides which task to execute next. 

To be specific, the priority queue model by Barabasi 
[1| consists of a single fixed-length queue, filled with tasks 
each of which is assigned a priority value drawn randomly 
when it enters into the queue. Every step the task with 
the highest priority is executed and is replaced by a new 
task with random priority value. Upon execution, the 
waiting time r, that is, how long the task has sat (waited) 
on the queue, is measured. Waiting time distribution 
P{t) of the Barabasi model has been shown to exhibit a 
power-law tail for large t as 

F(r)^r-", (1) 

with the exponent a = 1 0) Q i conforming to the be- 
haviors observed for the e-mail, library loan, and website 
visitation records 

0, i, [M [HI- Besides the human dy- 
namics, however, due to the extremal nature of its dy- 
namics the priority queue model would bear implications 
also to disparate problems in extremal dynamics, such as 
the Bak-Sneppen model for biological evolution [l^] and 
invasion percolation Q . 
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Barabasi model purposefully simplified many aspects 
of potential importance in realistic human dynamics, 
serving as a starting framework on which various detailed 
factors can be embedded [H, [13, [H, [H, [13 ■ One impor- 
tant factor that was not been accounted for is the human 
interaction. In the modern society, human engages in a 
large array of interactions with other individuals in var- 
ious modes. As a result, typical activity of a person is 
not an outcome of completely autonomous decisions, but 
of delicate compromises and balanced conflicts between 
often competing priorities. The impact of such a human 
interaction on the patterns of human dynamics, the wait- 
ing time distributions in particular, has been addressed 
recently by Oliveira and Vazquez (OV) [13]. They in- 
troduced a minimal model consisting of two interacting 
priority queues with interacting (I) and non-interacting 
(O) tasks. The human interaction is taken account for 
in a way that the /-task is executed only when both of 
the individual choose to execute them, that is, an AND- 
type protocol for the execution of /-task. Through this 
model they showed that the power-law waiting time dis- 
tribution still persists against the introduction of human 
interaction, but it has an effect that the power-law expo- 
nent a of P{t) can take numerable values other than 1 
depending on the queue length. 

Yet, the effect of human interactions for a system of 
more than two queues, or the queue network in general, 
has not been fully addressed. Such a question should be 
highly meaningful given the increasingly active engage- 
ment in various social networking of individuals forming 
complex network structures [Tsl [l9l. [20|. In this regard, 
here we study the dynamics of priority-queue network 
by generalizing the OV model. Specifically we focus on 
the scalability of the interaction protocols and the wait- 
ing time distribution under various model settings such 
as the number of queues (network size) and interaction 
topology. 

Priority- queue network — To construct a priority- 
queue network of N queues, we follow the OV model 
to divide the tasks into two classes; /- and O-tasks. A 
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queue has one 7-task for each neighbor in the network, in 
addition to an O-task. Thus a queue node i with degree 
ki (degree is the number of Hnks connected to the node) 
has a queue with fixed length Li = fc^ + 1 . We denote the 
/-task of the node i paired with the node j as lij, and 
the O-task of the node i as Oi. Given the network config- 
uration and the queue disciphne such as the interaction 
{e.g., AND- or OR- type) and update protocols {e.g., par- 
allel or sequential) , a priority- queue network is specified. 
Initially each task is given a priority value drawn from a 
uniform distribution in [0, 1). Then each step, each node 
chooses its highest priority task. The execution of the 
selected tasks is determined by the queue discipline: In 
the AND-type protocol, random sequential update case 
(the OV model), for example, we choose a random node, 
say i. If the highest priority task of i is Oi, then it is 
executed. If it is an /-task, say lij, it is executed only 
if Iji is also the highest priority task of the conjugate 
node j. Otherwise, node i executes Oi instead. The 
waiting times of the executed tasks are recorded, and the 
executed tasks are replaced with new tasks each with a 
random priority value in uniform [0, 1). N such updates 
constitute a Monte Carlo step (MCS), which is the time 
unit of waiting time measurement. 

OV model on networks — We first consider the general- 
ization of the OV model for N > 2 queues. We consider 
the model on two representative network configurations, 
the star graph and fully-connected network, for various 
TV. The resulting waiting time dynamics reveal an impor- 
tant phenomenon, the dynamic freezing due to priority 
conflict. The priority conflict occurs when a node i has 
lij as highest priority task, but the node j has another, 
say Ijk, as its highest priority, in conflict with each other. 

The star topology is less vulnerable to such a dynamic 
freezing since leaf-nodes can resolve it, primarily by up- 
dating priority of the O-task repeatedly. As a result, we 
have a power-law decaying P{t) (Fig. la-b). The power- 
law exponent a is found to be independent of the net- 
work size N; for the /-tasks ai w 2, and for the 0-tasks 
ao ~ 3 irrespective of being hub or leaf nodes. This re- 
sult is consistent with the OV model with L = 2, so on 
top of star graph, it behaves essentially the same as in 
the binary OV model. 

On the other hand, the dynamics is quite different 
in loopy networks such as the fully-connected networks, 
which are highly susceptible to conflicts that cannot be 
resolved readily. As a consequence, the number of exe- 
cuted /-tasks, r]{t), decays rapidly in time, either alge- 
braically for small N, or exponentially for large N > 10 
(Fig. Ic-d), and eventually the dynamics gets frozen, with 
the time scale decreasing with N. In real social networks, 
we have strong empirical evidences of high propensity 
of transitive triad relations (cyclic interactions between 
three individuals) f2Cil. 211 and clique (fully-connected 



subgraph) structure 20. |22|. so the AND-type interac- 
tion protocol would strongly suffer the dynamic freezing, 
rendering itself unrealistic towards realistic modeling of 
the network effects in human dynamics. 
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FIG. 1: (Color Online) (a— b) The waiting time distribution 
P{t) of the OV model on star networks for /-tasks (a) and 
for O-tasks (b), with various A'^ = 3, 4, 5, 20. Both P(r) decay 
with an asymptotic power law with A'^-independent exponent, 
ai ~ 1.9 (a) and ao ~ 2.8 (b), indicated respectively with 
dotted hnes. Deviations from the power law for large r for 
large A'' are due to the finite simulation time (10^° steps), (c— 
d) The double logarithmic (c) and semi-logarithmic (d) plots 
of the number of executed /-tasks, ri{t), for the OV model 
on fully-connected networks starting from random initial pri- 
ority assignments versus time. DiflFerent dotted line patterns 
are used for different network size A'^ (see legend). For small 
A', g{t) decays algebraically with A'^-dependent exponent (c) 
while it decays exponentially for large A'^ ^ 10 (d). Indicated 
slopes of dotted line are 1.5, 2.5, and 3.5, from right to left, 
drawn for the eye. Also shown is ri{t) for A' = 20 with the 
star graph topology (black solid) for comparison. 



Priority- queue network with OR-type protocol: The 
OR model — Not all human /-tasks should follow the 
AND-type protocol. As an alternative, an OR-type pro- 
tocol would be more reasonable for the tasks which re- 
quire simultaneous actions of two or more individuals but 
the action can be initiated primarily by either of them, 
such as the phone call conversation and the instant 
messaging [24'|. For such class of /-tasks, the potential 
priority conflict can be instantly overridden; we normally 
just answer the incoming phone call, for example. To 
model such situations, we introduce the priority-queue 
network with the OR-type interaction protocol by modi- 
fying the OV model as follows: a) Each step we choose a 
random node, say node i. b) If its highest priority task is 
an /-task, say lij , the two tasks lij and Iji are executed 
regardless of the priority value of Iji, If Oi is the high- 
est priority task, only that is executed, c) Priorities of 
all the executed tasks are randomly reassigned. We refer 
this model to as the OR model hereafter. 
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FIG. 2: (Color online) (a— b) The waiting time distribution 
P(r) of the OR model on star topology with TV = 3, 4, 5, 20. 
(a) P{r) for /-tasks decay as power laws asymptotically and 
the power-law exponent decreases as A'^ from qj ~ 3 for A'' = 3 
to ai ^ 1.5 for A'^ = 20. (b) P{t) for O-tasks show distinct be- 
haviors between the hub node (open symbols) and leaf nodes 
(full symbols). For the hub node, the power-law exponent 
varies from ao,hub ~ 1-5 for A'^ = 3 to a « 1 for A^ = 20. For 
the leaf nodes, -P(r) decays faster, with exponents ranging 
from ao,ieaf « 2 for A?^ = 3 to ao,ieaf ~ 1.5 for A'' = 20. 
For 7-task and leaf nodes' O-task, we shifted P{r) curve ver- 
tically to enhance visibility, (c— d) P(r) of the OR model 
on fully-connected topology, (c) P(r) for /-tasks decay with 
asymptotic powers, with exponents decreasing with N from 
a/ « 3 for AT = 3 to a/ « 2 for AT = 20. (d) P(r) for O- 
tasks follow the power-law decay with A'"- insensitive exponent 
ao ~ 1.5. Deviations for large r are due to finite simulation 
time. All quoted slopes are indicated with dotted lines drawn 
for the eye. 



For the OR model, P{t) still exhibit power-law tails 
for both star and fully-connected network topology, yet 
the power-law exponent a depends on the network size TV 
as well as the network topology in a diverse way. First, 
in the star topology, a decreases as TV increases: For 
/-tasks, it exhibits values from a/ « 3 for TV = 3 to 
ai « 1.5 for TV = 20 (Fig. 2a); For O-tasks, the expo- 
nent exhibits distinct values for the hub and leaf nodes, 
changing from ao.hub ~ 1.5 for TV = 3 to ao.hub ~ 1.0 
for TV = 20 for the hub node, whereas for the leaf nodes 
it changes from ao.ieaf ~ 2.0 for TV = 3 to ao.ieaf ~ 1-5 
for TV = 20 (Fig. 2b). Moreover, for the hub- node, the 
mean waiting time {T)o,hub of O-tasks diverges with the 
exponent a < 2, similarly to the Poisson queue placed 
on the hubs in scale-free networks [l^. For other tasks, 
P(r) with a < 2 for large TV is accompanies by the peak 
at r = 1, rendering the average waiting time finite. 



FIG. 3: (Color online) The waiting time distribution P(r) 
of the OR model with parallel updates, (a— b) P(r) of the 
parallel OR model on star topology for the /-tasks (a) and the 
O-tasks (b). For /-tasks, the power-law decay of P(r) exhibits 
a weak TV-dependence, as the exponent varies from ai ~ 2.2 
for TV = 3 to ai « 1.7 for TV = 20. For O-tasks, hub and leaf 
nodes exhibit different TV-insensitive power laws with aj^hub ~ 
1.5 for the hub and aj^i^af ~ 2.5 for the leaf nodes, (c— d) 
P(r) of the parallel OR model on fully-connected topology. 
For /-tasks the power-law exponent is found to be insensitive 
to TV as a/ ~ 2 (c). For O-tasks, however, the asymptotic 
power-law exponent increases with TV, from ao ~ 2 for TV = 3 
to ao ~ 3 for TV = 20. In (d), the P{t) curves for different 
TV are shifted vertically to enhance comprehensibility. All 
quoted slopes are indicated with dotted lines drawn for the 
eye. 



In the fully-connected topology, the power-law expo- 
nent a weakly depends on TV. For the /-tasks, it de- 
creases with TV from a/ « 3 for TV = 3 to a/ « 2 for 
TV = 20 (Fig. 2c). For the O-tasks, on the other hand, 
ao is rather stable against TV as ao ~ 1.5 (Fig. 2d). 
This result implies that on the fully-connected networks, 
/-tasks are executed with finite mean waiting times while 
O-tasks on average have to wait on the queue infinitely 
long to be executed. Taken together, the OR model on 
networks implicates the importance of not only the over- 
all network structure but also individual node's topolog- 
ical position on determining the dynamics of networking 
priority-queue nodes. 

Priority- queue network with parallel updates — Update 
rule in discrete time dynamic models has been known to 
affect the dynamics considerably [2^ . Thus it is informa- 
tive to study the effect of update rule in priority-queue 
models. To this end, we consider the parallel update rule 
by modifying the OR model as follows: a) Each step, 
each node chooses its highest priority task, b) We sort 
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all chosen tasks by the priority values, and execute them 
in order of priority, while each node can execute at most 
one task each step. That is, if the priority of /y is higher 
than that of Ijk, than the node j executes the task Iji 
upon request from i before Ijk, which subsequently can- 
not be executed in this step, c) All the executed tasks 
are assigned new random priorities, completing a MCS. 
Priority conflict may occur at the step a), but it is not as 
problematic as in the AND-type protocol case even in the 
fully-connected topology, due to the partial resolution in 
b). We found, however, that it is strong enough to affect 
the waiting dynamics of the priority-queue network: it 
can reshape P{t) in a significant way, because the tasks 
in the tail (r 7^ 1) are particularly strongly affected. 

The waiting time distribution P(t) of the OR model 
with parallel update unanimously exhibits a probabil- 
ity weight strongly concentrated at t = 1 followed by 
an asymptotic algebraic tail (Fig. 3). On the star graph 
topology, the decay exponent a shows less variation com- 
pared to the sequential update case: It changes from 
ai « 2.2 for TV = 3 to a/ ^ 1.7 for TV = 20 for the 
/-tasks (Fig. 3a). For the O-tasks, it is insensitive to N, 
yet exhibits different values for the hub and leaf nodes, 
as ao,hub ~ 1-5 and ao,ieaf « 2.5, respectively (Fig. 3b), 
decaying faster than the sequential case. For the fully- 
connected networks, the power-law decay exponent for 
the /-tasks is a/ « 2, insensitive to N (Fig. 3c), while 
for the O-tasks, it even increases with N, from ao ~ 2 
for TV = 3 to ao w 3 for TV = 20 (Fig. 3d). The presence 
of strong peak at r = 1 renders the mean waiting time of 
all tasks finite, which is due to the partial resolution rule 
we implemented in the model. These results demonstrate 
clearly that the update rule in the priority queue network 
strongly affect the overall dynamics in a nontrivial way. 

Summary and Discussion — In this paper, we have gen- 
eralized the binary interacting priority queue model, the 
OV model, into a queue network, showing that the OV 
model is not easily generalizable onto loopy networks as 
the dynamics gets frozen due to substantial priority con- 



flicts. We then introduced a modified model, the OR 
model, which can be put on top of any network topology. 
It is shown that OR model exhibits power-law decaying 
waiting time distribution P(t), Eq. (1), yet with diverse 
values of the exponent a depending on the global network 
topology, local position of queue nodes on the network, 
as well as the queue discipline such as the update rule. 

The fundamental factor driving the diverse behaviors 
of the interacting priority queue models is the existence 
of priority conflicts. Different global network topology 
and distinct network position of nodes impose different 
degree of conflicts and the resolution thereof. In the per- 
spective of human dynamics modeling, its introduction 
seems reasonable, for such conflicts are more of a rule 
than an exception in daily decision making of modern hu- 
man life. A crucial question remains at this point. The 
dependence of waiting time dynamics of a priority queue 
node (a human) on the global or local network topology 
that we found in this work has never been gauged from 
the empirical data yet. Candidate datasets for this end 
would be the mobile phone data fgsj or the instant mes- 
saging data .24 ], for which power-law- like waiting time 
distributions have been reported. The waiting time t 
in this work is measured for each task, that is for each 
link, differently from what have been previously mea- 
sured from data in Refs. [23|, [l^], so a direct comparison 
cannot be made. Appropriate measurements with these 
datasets would reveal the relevance of the priority-queue 
network models studied in this work and the role of in- 
teractions in human dynamics in general. 
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